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ABSTRACT 

The kinetic model for irreversible novolac type for 
phenol formaldehyde polymerization is presented and equations 
governing the molecular weight distribution MKT) of the polymer 
formed for both batch and homogeneous continuous flow stirred 
tank reactors (HCSTR) have been derived. The set of non-linear 
differential equations for batch reactor is solved using Runge— 
Kutta method of order 4 and the set of non-linear algebraic 
equations for HCSTRs is solved with the help of Brown's algoritl-irn 
which is found to be more efficient than Newton's method, A 
sensitivity analysis of different reaction parameters hcive been 
carried out and the reactivity of the para position was found 
to be an important factor affecting the 141)®, The results for 
HCSTR have also been compared xvith those for batch reactors and 
the former is found to give polymer of lower average molecular 
x\reights having higher polydispersity indexi A model for the 
reversible case is also presented. 



CHAPTER 1 


Introdu ct ion i 

There are several ■ classes of condensation polymerization 

systems where Flory’s equal reactivity hypothesis does not hold 
2 

good • Formation of Novolac type phenol formaldehyde polymer 
is one of the commercially important systems. 

It is known that formaldehyde in x^ater exists as 

methylene glycol (OH - CH2“0H) and therefore exhibits a functional! 

of two* Phenol on the other hand has three reactive sites^ two 

ortho positions and one para positions and experimental studies 

^'3 

have shovm that they have different reactivities. Also the 
studies on gelations , have revealed that the isoliTiier molecules 
grovj mainly by the reactions of the end groups* Drum and LeBlanc 
have explained this phenomena through molecular shielding assigning 
lovjer reactivities to internal sites*. 

For polymerization systems where equal reactivity 
hypothesis holds, the expression for the reaction rate can be 
written in terms of reaction between functional groups* But in 
the case of novolac type polymerization one must write mass 
balance equations of individual polymer molecules of a given 
chain length in terms of reaction of these site;^and the resulting 
system of differential equations (including the balance equations 
of different sites) must be solved n^ameri cal ly for getting 


conversion and MWD 



CHAPTER 2 


Kinetic M.odel ; 

In the reaction mass, the polymer chains are assmed 
to consist of all those molecules X'/hich have chain lengths greater 
than 2. There are txTO kinds of molecular species 



x/here the s;i'mbol means that the bond 'formation has occured 

i 

either at para or ortho position. The p and Q differ fron each 
other only by substituted -CH 2 OH, It is assiomed that there is 
only one bound CH 2 OH per chain which can be anyvjhere on it. These 
molecules contain external ortho (O ^p) and papa loosition (p ) 
which are present at chain ends and internal ortho and para 
positions (0^ and ,p^) situated uithin the chain. All these 
positions have been shovm to differ in their reactivities. The 
examination of chain ends reveals that the following two molecular 
structures can give rise to external ortho positions 
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These are distinguished and denoted as 0^ and 0^* The total number 
of external ortho positions is the sura of these two. 

At a given time in the reaction mass, there can be 
free formaldehyde molecules and bound CH 2 OH groups and the 
reactive sites can interact with both of these. The different 
reactions have been siaramarized in table 2.1, These reactions haive 
been vrritten to satisfy the site balance. As an examx:>le, reaction 
(a) table 2J. stands for the fact that on reaction y 0^/ is assumed 
to be not available for further reaction v/hich becomes on 
further reaction. This fact has been rexsresented by writing- 
-O and +0. in the right hand side. The reactivity of these 
reactions is assumed to be ccmpletely determined by the reactive 
sites involved in them. In reactions (a) to (e) of the table, 
the rate constants have been multiiDlied by factor 2 because in 
the fon'naldehyde molecules, there are two (OH) groups where the 
reaction can occur vmth equal likelihood. 

The schematic reactions shown in Table 2 JL do not include 
these involving phenol. When phenol reacts with formaldehyde, it 
forms the following two species: 


OH OH 



CH2OH 
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Reactions of and are also not covered by schematic equations 
of Table 2*land they are all listed in TaJole 2»2, In reaction (a) 
of table 2.2 there are two positions on iDhenol (P) and two 
positions on F and therefore the reactivity has been shovjn to be 
4k^ and so on, A mass balance on these sioecies for batch reactors 
gives 


dt 


dt 


coil 


[QJ + 2 [f1 + [CHjOh] . + Toy 


+2^2 [Oj +2k3 t P^l +2k^ [ +4(k3+ky [PJ 


[Ojl +L«2 °h) I 2k^ [P] 


(4) 




= thl 


^1 L°41 + h [°ll ■^’^3 [Pi] 


’"iW •"’^4 [Pel [Pi 


-[Qyx 


f>l 


kp [0^1 


4 r Pe] 


f CH2 OHI 


^4 fPe] [O^} 


(5) 


r n 1 


-cSt t“ll I [°4l ->^2 LOil 41^4 f Pel! 't. fe'l - >^2 [°li 

L _ __ I_ 

■‘■^4 C Pel 1 ■''2 [ f] I [ 0 ^ J - ^=2 L°ll '*'^4 [Pel 


4 fcHj Oh 3 jK^O^l -1^2 [Oj + k^ [p^] 


(6) 
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at- - [°el I l°l] W + ' -■ 

i- •' 

I [Qi'i + [Oil + 2 [f] + [cHjOh] 


2 °H] 1 -=^3 I Pi ] 


(7) 


aJPel 

dt 


L°lll ""l L°4l [°ll +(41 ^i+K 4) I>] +K3tPi]J + [Qilx 

! 2ic, [P] -k, [Oj -k^ (p^ jl -2 Tf] ' k^ [0^1 +k^ fpj 

t- 


+ [cH20Hl[2kj [p] -k^ lp^lj 


( 8 ) 


-Tt— = b] - 2 (ki 4 * 4 ) tQilW -L°ll t°il +>^1 C°al 

■^^^2 L°il "^^^3 1 ^ 1 ] ■^’" 4 [Pg 1 \ -> 2 (]c^+]c^) [Q^"] -( 3 k^+]c^);>^ 


[Qil [>il “Oc^+k^) [ch^oh] 


(9) 


^T“- = W -(Skj^+k^) [Qilf Pj -tQ{l 


>^1 C°i] 


tk^LOil -•■k3 [Pj +k^!:PelJ -4>^1 CQi]"-(3V’^4> C°l!, 


l°i] [Qll [cHjOh] 


(10) 


d L f] 




(2kj+k^) [P] +k^ [ 0^3 +kj^ [0^3 +k 2 [Oj^] ■HC 3 


[ Ppl + k^ [p^- 


( 11 ) 


d 


(2k^+k4) C P] + Coy + 2 Cf) + CCH 2 OHI (12) 
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Modelling of Molecular Weight Distribution : 

In previous studies^ ±t leas shox^n that the concentration 
of bound CH 2 OH in the reaction mass is small. In viev7 of tliis^ it 
is assumed that- a polymer molecule in the reaction mass either has 
one bound GH 2 OH or has none, .In view of this, the grovjth of 
polymer molecules can be represented by following reactions: 


and 


P -f- F 
n 


n 


Q + P - 
n m 


^n 


' n+m 


Q. 


'n+m 


n,m =2,3 

n,m =2,3 ,,, 


(a) 

.(b) (13) 

(c) 


and 


The reaction batv;een and can occur in two ways: 

(a) reactive sites on react with (-CH 2 OH) on 

(b) reactive sites on react with (-CH 2 OH) on 

If O^Q^, ®i^n' ^i^n ^e^n total number of 

sites on all in the reaction mass and *^i^' ^i^' 

on all in the system, the reaction rate, betvjeen 

and can be written as 


mn 


= “n ’"s % “n + C 

+ 3n * '"2°1 + >"3 n 3a + 3a> L 2n] 

m, n =3,4 ( 14) 


The a-bovB equation must be modified vjhen either m or n is equal to 2 
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This is because polyraer molecules of chin length 2 do not have 
internal sites and the above equation is modified to 


°e Sn + >^2 ^ [ °2J 

+ (k^ + kjO^ Qj + op [c^] m = 3;4 (15) 

The reaction between Q and P , hov/ever, xirould not obey 

n m 

eqxiations (13) and (14) because the bound CH^OH of Q alone can 

^ in 

react vrith reactive sites of P^. Thus, the corresponding rate of 
reaction, can be xirritten as 

r‘ = (k. O P + k,0' P + k„0. P + k-, p. P 
mn l em lem 2im 3^im 

+ k,p P ) Tq 1 n = 2,3 
4-^e m t“ni 

m = 3,4 (16) 

When m = 2, the above equation must be modified to 

°G ^2 + >^1 °i ^2 ’^4Po P2> 

To find Q^, it is noted that 

= f Q 1 (Average fraction of sites on Q„ that 
G n ^ n n 

are 0^ ) (18) 

G 

The fraction of sites on is assumed to be the same, on an 
average, as that of the overall average in the reaction mass. 

In other words, 

O Q -Q: fo 1 (overall fraction of 0 sites on all 
e n >'■ n.' e 
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polymer molecules) 


= fO. i’ 


\oJ 




t. n 




I 


(19) 


Following relations are written for other reactive sites as 

t°i ! L°nJ 


O* Q 
G n 


O.Q 

in 


p . Q == 
*^1 n 


P Q 
•^e n 


SUM 


Ihi l^nl 


SUM 

Lpi ) 

SUM 


[°n.l 


I Po l L^n! 


.» Lh 
SUM 


where 


SUM = 


(Z. 


n 


+ -ii-' Q„) 
2 ^ 


(a) 

(b) 

(c) 

(d) (20) 

(21) 


Similar relations are valid and can be written for 0^ 

etc. In terms of these^ equations (14) to (17) can be easily 
written as 

L L Q-’i 


2 CN 


mn 


ma 


SUM 

CN 

sum'"' 


n ’ 


[y [°2l suf^ lAlM 
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q 


mn 


2n 


Hr' '>mlW 


CM 


SU14 


bal [ O 


(a) 

(d) 


(22) 


whoro 


CN - k^[ O^j + {oil + 1^2\6±1 ^"3 l^D. ^4 L^ei 


and 


i-°el ^1 t-^el ^4 


(b) (23) 


Now one is in position to cariry out the mass balance 

on P and Q which is done as follows. For the mass balance of 
n n 

it is noted that it is formed when reacts vrith formaldehyde 

and Q reacts with Q _ v^hero m could have any value between 

m n— m: 

1 and (n-1), is destroyed v;hon it reacts with phenol or P^ 

and where m =2^3,,.,, The balance for a batch reactor yields 


= -(2kj4.kp i_p1Lq„) - ^CN [Q^i + -gn + [Pii 

( m=2 

+ ap) + (2k^ [q£] + 

" 4*' 2 [hi [VJ - S^[°2l [V2> ^ 

i=2 

Lhl [h-3 Ph Pn-fj + si- CVll Pi] 

+ [ Q •] n > 4 


( 24 ) 
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The relation above is valid for ny^4. and the balances for Q„ and 


are 


<3 L 


dt““ = COil +3(Krk4) CQil tap CQif+ ffzlhl 


-(2k,+ic^) (p-j LP -CN LQ2l(?Lag + LQj+bp) 


tP 2 -i (2 k ^[ Q |3 + L^il ■*■ ^'^4 11%^ ^ 


(25) 


and 


d CQsl 


(2ICi+k^) CQ3-3LP] -CN rQ3T - si- LQ3] ( + top 


dt 


Z 


+ 


2 CN 


+ (k^+k^JLQ^] [“al+^h *^°2^ 


L^i ^ ^ ■ SUM 

+ s^"- !;a2l <■ toi'i + CQp ) CQp +)^1 LOil +k 

LQil ) fQ3] 


4 

- (26) 


On summation of equations (24) to (26)/ one obtains 

d i:-ch20h] d |I[; op 


CXk. 

TT 


OL 


dt 


•f 


dt 

2 £f3 


SUM 


- = -(2k^+k4) L”P ^V^Qnl 

-ns*- 

(0^2 LPj + CMpP^l ) + (kj+k^) t_Qj) 


yt-i 


+ Okj+k^) Co^ [Qp + 2k cop 


(27) 


A smilar ma ss balance can be made on It is noted 


that P is formed when P . reacts \fith Q- or Ql and P . reacts "w-, J 

n n—l 1 1 n^i 
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It is dGstroyed \7hen it reacts with formaldehyde and any 
The balance gives 


<sCPn3 


n-3 


dt 


SUM 


Lti-il < Ph ^ PM S si--LvM 

i=2 


sS- Lhl IVal L^SLvM - PM 


CN 

SUM 


Pnl (ZZltg + [Pil + Pi] > 

2 


(28) 


The above equation is valid only for n 5 and therefore balances 
for P 2 / P 3 and are written as 


a Up,] 


dt 


( 2 k^+k^) [p] + j 5 {.]) - 


" “ Ld 


SUM 


CM 

SUM 


LP 2 I < ^ ^ 

2 


(29) 


1113 . 

dt 


cm 


=- 1^21 * [Pil )+(2ki+k4) ilPJ CQj) - “1^ LPaiiF 


SUM 


CN 

SUM 


Cp 31 ( Z.0^1 + Top + iPll) 


d 


i:p 4 ) 


CN 


[Psl^ 1P{1 ^ ^2^ ’*■ 


dt 


SUM 


X P3] - -4® PMP] PM ' 2:t.cwi + cojiH-coi.,) 

^ (31) 


CN 


SUM 



12 


To calculate the molecular weight distribution/ eqns. 

(4) to (12) and (27) must be solved simultaneously and substituted 
in eqns.- (24) to (28)/ to obtain fP I and \ Q \ as a function 
of time. These equations are written in dimensionless form 
using following variables* 



(a) 


^2 - IFio 

(b) 

y3 = ITo^l / &■] „ 

(c) 

^4 “ / h] o 

(d) 

Ys = CPe'] / I>1 o 

(e) 

Xs = LQi] /Cp] o 

(f) 

Y-, = CQ{V [Pi o 

(g) 

Ys = Lp^/C P]o 

(hi) 

Yg = / r P] o 

(dt) 

y^o= / tP^ 

(j) 

yii= C PaT / ilpJo 

(lO 

^ = h O' ! 

(1) 

= tPil/LPlo*” 

(ra ) 
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^Q. = 

1 

CQi^i / 

(n) 

II 

(k2/k^) 

(o) 

R2 = 

(ka/ki ) 

(p) 

^^3 = 

(k^/lc^) 

(q) 


Ha mogeneous Continous Flow Stirred Tanlc Reactor f 

The rate equations derived for batch reactor can be 
incorporated in the steady state HCSTR simulation as follows: 
Rate of input + Rate of generation due to reaction 
= Rate of out put 


For P and Q 
n n 


or 


o'"'" (rp ) = V L'Pnl 

n 

fp 1 - ip 1 + r = o 

L n J o V- n-* dt 

[°nl o - t?n1 = ° 


where r = rate of formation of P^ 
n 

V = veluiwe flow rate V*|urv\c^/time 

V = volume of reactor volume 
Klc input concentration mole/volume 

C^l ,|p^l = out put concentration mole/ volume 
T = residence time 

Similar equations can be written for ail the sites> 


t 


P 


F 
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Q. 


, Qi and 


Thus vre get a system of nonlinear 


algebraic equations. Since equation for -CH^OH is a dependent 
one, it is not incorporated. We get -CH^OH frcm the relation 


C- CH2OHJ = 

n=2 





fract 


on distribution MFDP^ and MFDQ^ of 


The mole 

species P. and Q. resepctively have been defined as 


MJDP. + — .rr* 
1 ^ 




iitpd + 1>’ 
1=2 1 




MFDQ^ = 


rql 


' 2 - l Qjl + [oj; ] 

i=l 


r 1 / 2,, , , 


In the above equatij^n for n=l in MEDQ^, is equal to ( + Q| ) 

has been used in the calculation. 

In the weig'it fraction distributions, WFDP^ and WFDQj^, 
as opposed to this, the total weight of the reaction mass at the 
time of interest has been made as the basis as follows; 


WFDP . = 

1 & 


WP. 

1 


-2^ I/jp. +2^- WQ. + WP + WQf 
l=2 1 i=i 1 1 


and 


WEDQ. = 


WQj_ 


WP. + -2- WQ. + WP + WQ* 
i=2 ^ i=l ^ ^ 


where WP., WQ., WP and WQ* refer to the weight of species P.,Q.,P 

" ■ i ’ JL, ' X , .X, -i. ^ 


and Q^ respectively, 
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TABLE 2.1 

Reactions of O / O', p. , O. and p 

, , ... „ _X _^ 0 


(1) with formaldehyde 


21c, 


04 + p 




0^ + F 


Pj_ -!- F 


Pe + F 








2k. 


4 



^ (-CH 2 OH) + - ( 

(-CH2OH) + p^»pg 

“> (-CH 2 OH) 

■> (-CH 2 OH) 

^ (-CH^OH) + 0.~0 

^ J. 0 


(2) With bound formaldehyde (-CH^^OH) 


•^4 


O' + (-ai,,OH) - — 1 - .- ^ o.~( 

e ii 1 

ki 

°e (-CH2OH) p^-p^ 

^2 

Oj^ + (-CH2OH) — — — consumption 

^3 

Pjj_ + (-CH2OH) — ' >■ consumption 

Pg + (-CH2OH) 


°i - °e 


(a) 


(b) 


(c) 

(d) 

(e) 


(f ) 

(g) 

(h) 

(i) 

(j) 
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TABLE 2 .2 


Reactions of and Q| 


(a) Reactions of P with F and (-CH 2 OH) 


P+(-CH20H) 


Pe-^ Oe 


P+(~CH20H) 


2 O' 
e 


(b) Reactions of P v/ith and 


P + Qj[ 


p H- 20' +0 

•^e e e 


P + 


4 O' 
e 


P + Q. 


2p + 20^ 

■^e e 


P + Q. 


O + P +20' 
e e e 


(c) Reactions of and Q| with reactive sites 


O' + Q* 
e e 


o. + 20' - O' 

i 00 




p + O +0.-0' 
■^e e 1 e 


Oe + Qi 


Pe + V Pi~Pe 


^e + Q1 


20^ + Pj_*"PQ 


(1) 
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p . + Q-i 
1 


Pi 


0. + Qi 


O . + Q« 
1 1 


Po + Qi 


Pe + 




k. 




O + p 

e 

(m) 


(n) 

20- 

(q) 

o p 

(r) 

°1 •" °e+ Pe-°e 

(s) 

0 l+ 20 ^ - °e 




CHAPTER 3 


numerical solutions 

1 • Etunge Xutta Method of Order 4 ;' 

The system of nonlinear differential equations for the 
batch reactor are solved in DEC 1090 computer using Runge Kutta 
Method of order 4, After a fei^'r trial/ the stable solution was 
found for x = 10 , The algorithni of the technique is 
as follows^: 

For the equation 


yi 

= h 

(x. /X„/X_/ , 

1 A O 

*. *./X^) 

^2 

= ^2 

( x^/X2/X3/ 

n 


= ^n 

(x^/X2/X3/ , 



with y (Xq) = Yn v;here the bar ‘indicates array gnerate 


0 


apjproximations y^ to y ( Xq+ n ^ x ) for ^ x fixed and for n 
n=0/ 1/2/,.., using recursion formula 


n+l 


^n ^ \ +2}C2+2k3+k4) 


where 


k. 


X f (x ) 
n 

X f (x^ 4 - ) 

X f (X^ + } 

Bb 4 


‘4 


X f (x + 4X) 

■ n 
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2 , Brown's Method i 

The following system of nonlinear algebraic equations 
is considered 

f^ (x) = f^ (x^/X^/ . . ./-Xjj) = 0 

±2 (x) = ±2 (x^,X2/....x^) = 0 (3,1) 


‘N ° 


xdiich can be written in vector form as f (x) = 0, 

It is assumed that f^^ are continuously differentiable. In 
Ne\^on Raphson technique the functions in equations 3,1 are 
expanded about a point x^ which is close to the solution x* 
yielding 

f (x) = f (x^) + j (x^), (x - x’^) + H.O.T. (3.2) 

where the higher order terms (H.O.T.) are negligible and only 
the first order tenns are included. Since x^ are close to x* 
one therefore has 

f (x) =0 or 0 - f (x^) + j (x^) , (x-x^) 

Solving for x gives 

= X = x^ - (x^). f(x^) (3.3) 

The solution pro<2edure would therefore be the following: 

Begining with a starting guess, x one solves equations 3.3 
iteratively for n= 0,1,2,,,., This can be done only when J (x) 
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is non— singular in a suitalDle neighbourhood of x*. a derivative 
free analogue of Nex^^ton's metbd often called the discrete Newton's 
method is obtained fran equation 3,3 by replacing the ith row and 
jth column entry^ of the Jacobian matrijc (i.e, ,r— = ) by the first 
difference quotient aprjroximation 


f.(x^+h^.) - f. (x^) 
i 1 2. ^ ' 


n 


(3.4) 


Here e^ denotes the jth unit factor and scalar h^ is normally 
chosen such that h^ = 0 (11 f (x^) 11 ) . 


Brown technique calls the functions of equation 3,1 
one at a time so that information obtained from v/orking with f ^ 
can be incorporated in while v/orking v/ith etc, A successive 
sxfDstitution scheme is used rather than the simultanoous 
treatment of the f., as done in Newton's method. The former 
method is derivative free and its second order convergence has 
been proved. The method consists of applying the following steps; 
Step 1 ; 

Let x^ denote an apj/roximation to the solution x* of 
3,1, The first function is expanded in an approximate Taylor 
series expansion about the point x^. By approjciraate it is meant 
the expansion in which the actual (analytic) partial derivatives 
are replaced by first difference quotient approscimation. 
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f, (x)= f, (x^) + f 


' 1 ^ ' ^ ^Ix 

.n X . . ^ ,-^n 


x^;h 


rn, 


(3.5) 


Here ^ defined to be the first difference quotient 

x-j/ 


f2(^) 


approximation 3,4 with i=l. If x^ is close enough to x*^ 

O and equation 3,5 is equated to zero and solved for 

that variable, say x^^, whose corresxoonding approximate partial 

derivative, f. (x^) is the largest in absolute value. This 
ix^,n 

gives 






N' 


where ^ix^yhi given in 3,4 and f^ = . 


The constants f^^ f^x ‘h saved (stored) in the ccxnputer 

implementation of the algorithm for future use , By choosing the 
approximate partial derivatives of the largest value to divide 
by a partial pivoting effect is achieved, sirailar to what is often 
done when using the Gaussian elimination process for solving 
linear systems. This enhances the numerical stability of the 
method. It is observed frean equation 3,6 that x^ is a linear 
function of the N-1 variables x^,X 2 /*.. .x^ ^ and for the purpose 
of clarity, the left hand side of 3,6 is redefined as 


-n 


L^(Xi,X 2 / ....x^_^l) and define (x^^x^„#,, •^»1> 
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Step 2 ; 

Now a new function of N~1 variables ^ 

is redefined which is related to the second function £2 of 
equation 3,1 as 


Let g 2 be defined as = £2 (^ 1 ' '^2 

expanded in an approximate Taylor series exxoansion about the xsoint 

(x^, ,.,x^_^)/ lineari2ed (ignoring higher order terms) and solved 

for that variable^, say x^^ whose corresponding axaxjroximate 

laartial derivative g^ , is largest in mangnitude, 

^^N-l>h. 


X, 


N-l 


n 

"n-1 


N-2 


j=l 


(g 


n 


X . ;h'^^2x. 
3 


N-l 


;h>- 


;h 


(3.8) 


Here the approximate partial derivative g^^ ^ is given by 


n 


/ n n n , . n n ■ n ^ n 

92(==i q.i'q 


^2x .;h 


n 


(3.9) 


where h^ is a small positive number, Fran equation 3,8 it is 
noted that x^ ^ is a linear function of the remaining N—2 
variables and it is denoted by ^(x^z ,,,/X^__ 2 )* Again the 

ratios g2x^.;h^^2x^__^;h , j=l,..,,N-2 and gJ/S^x^ 
for future use. 


are stored 
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Step 


. 


One defines 


(^- 1 1 * t ' 


) f ^ (x. / « * O /^TVT 1 


’3 "■3^“l'-“-^^N-2rN-l'-"^N' 


(3.10) 


with the arguraent of being (x^, , , ,,x^_ 2 ) and the argument of 


being 


^1'^2'***'Xn_2' ^n-1 •••^^N-2^ 


(3.11) 


One carries out 1) approximate Taylor expansion of the 

function about followed by 2) linearization of the 

resulting expansion follovred by 3) equating to zero and solving 

for one variable^ say 3c^__2 (whose corresponding approximate partial 

derivative g_ u is largest in magnitude) as a linear 

■^^N-2'“ 

combination 2 / of ihe now remaing N-3 variables. It is 
continued in this fashion, eliminating one variable for each 


equation treated. Every time a nev7 linear expression, L.^ , , 

N*-K 

one of the variables, say x_. , , is obtained in terms of the 

N—K 


for 


remaining N-k-1 variables, x^,X 2 ,...Xj^ 2'^N k 1* This linear 

expression is used whereever x^^ had aiz^ioeared in the previously 
defined linear expressions , ,,, L„ , It may be noted 

that each step in the alogrithm adds one more linear expression 
to a linear system. During the k+1 st step of algorithin, it is 
necessary to evaluate 9^^ 2 . ^]c+i' various arguments. 

The values of the last k components of the argument of 3-re 

obtained by back s\ibstitution in the linear sy stem L. ,, L„ i, , . . WE. , . 

INI JN**Jv*r 
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i^hich has already been found* The points which are back 
B-ubstituted consist of the point (x^^,..*,x5 , ) = xJJ together 

JL N"»jC N*^JC 

■;ith the points + h^e^, j=l,...,N-k, v/here e^. denoted the 

jth unl't vector* The argijraents are required to determine the 
[uantitities and ^ j=l^ .,.,N-]c, needed for the 

ilimination of the k+1 st variable; say x^^ by the basic 
processes of expansion, linearization and solution of the resulting 
;xf)ression. This x^tocess results for each k, in the k+1 st 
•ariable, say x , > being exxoressed as linear corabination' L.^ , 
if the remaining N-k~l variables, 
tep N: 

I33.> the last step 


= f._ (x- . ,,L-.) where L‘s are obtained by bacJc 
15 1 iM J 

ubstitution in the N-1 rowed trangular linear system which now 
as the form 

i“l 




j=i 


^%-i+l,x ^ ^ ^N-1+1'^ 


Cl+l,x.;h 


(3.12) 


arith g^= f^ and = x^) so that g^ is just a fionction of the 
Ingle variable x^. Now expanding, linearizing and solving for 
, one obtains 
X, 


L 


^1= - 9>S,x,,h 

lis point x^ thus obtained is used as the next axapraxiraation. 
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X i”T“ J_ 

, two the first component x*, of the solution vector x*. 
x^ is renamed as and the system in equation 3,12 is solved 
to get improved approximations of other components of x*. 

The successive substitution nature of the algorithm 
allows the most recent information available to be used in the 
construction of the next function argument similar to what is 
done in the Gauss— Seidel process for linear and nonlinear 
systems of equations. A/ simplified flow chart of the Brown's 
method is shov/n in figure 3,1 



formulate g" in 
terms of Xj^'s .Lj^'s 



Guess a solution 


Calculate fDarfiol 
, derivatives 


Maximum partial 
derivative 


/Partial pivoting based on the 
\ above to choose thevoriabi<^ 


Solve the equation having 
one unknown and 


Bock substitute to get 
V n*>t 1/ ^ o 










CHAPTER 4 


RESULTS A]MD DISCUSSION 

The mass balance equations for both batch as well as 
HCSTR are solved numerically in a DEC-1090 computer to examine 
the effects of various system iDarameters on molecular weights 
and conversions. 

In figure 4,1 and 4,2 the mole fraction distributions 
(MED) of and have been plotted for different residence 
times of polymerization in HCSTR, On the same graph, the 
corresponding results for batch reactors have been given for 
comparison. As the time of polymerization increases, the 
distributions are foxind to become broader. For any given x , 
the conversion of phenol is always lower and (= + LQ^ f ) 

is higher for the HCSTR, In view of lower conversion of phenol, 
the concentration of higher homologs are naturally lower than those 
obtained in batch reactors. In figure 4,3 and 4,4, the 
corresponding weight fraction distributions (MWD) for HCSTRs as 
well as batch reactors are given. The distribution for batch 
reactors in figure 4,3 undergoes a maximum at n=2 for residence 
times beyond x=0,3, whereas no such behaviour is seen for HCSTRs, 
The distribution curves for p. as well as Q. for HCSTR's are 
always found to be lower than those for batch reactors. In 
figure 4,5 the number average^ and weight average, M^, 

molecular weights of polymer formed (including both and 
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in the reaction mass for batch as v^eil as HCSTRs has been plotted 

as a function of time, as expected, M and M for the latter are 

n w 

always found to be lower; however the polydispersity inde^c, , 
for HCSTR is larger. 

In figure 4,6 and 4,7, the effect of R^ on the mole 
fraction distribution of and at dimensionless residence 
time x=0,9 has been examined. In figure 4,6 for a given R^/ the 
distribution of falls monotonically and beyond n=10, it -is 

almost zero and as R^ is increased, it becomes slightly broader. 
As can be seen from the figure the effect of R^ on the MWD of the 
polymer from batch reactors, in contrast to this, is found to be 
much more. In figure 4,7 the effect of R 2 on the mole fraction 
distribution of has been shown, LQ^ l + LQ|1 ) for 

HCSTR is found to be much lQ3|iger than for batch reactor, whereas 
the remainder of the distribution curve is fomd to be negligibly 
affected. 

In figure 4,8 and 4,9, the effect of R^ on and M 
versus time has been examined. As R^ increases, the curves 
rise upwards; however for a given R^ they are always lower than 
those for batch reactors. It can be noted that for large 
the reaction at para positions occur xirith higher likelihood, this 
way the preferred molecular structiore being OFT — CH2~^^/‘0H, 
This leads to the reduction in the polydispersity and gives 
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considerably lovKsred M , For the values of studied for HCSTR,' 
in figure 4,9^ this is not found to occur even though the cuiTvres 
become close to each other for R^ beyond 5, 

In figure 4,10 and 4,11, the effect of CpJ on 

the mole fraction distribution of. and have been exai'ained. 

In figure 4,10, the conversion of phenol (for x=0,20) is always 
foTond to be smaller and as a result of which the distribution 
curves appear steeper than those for batch reactors. Similarly 

(= HCSTR is larger. On 

calculation of the polydispersity index, , of the polymer 
found in HCSTR is always larger. In table 4,1 and 4,2 the effect 
of R^ has been examined on the distributions of Pj^ and and 
results compared with those for batch reactors. Similar results 
are obtained when R 2 is varied and the MWD are foiond to be 
insensitive to these parametrs. 


in figure 4,11, 
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FIGURE 4,1-- Effect of Residence Time on MFD of P 




PIGUR 4.2 ~ Effect of Residence Time on MPD of 
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CHAPTER 5 


MODEL WITH REVERSIBLE KINETICS 


Earlier it was ass-umed that the condensation product 
water is removed continuously and hence reverse reactions donot 
take place. When this is not the case the kinetic model must be 
modified, which is done as follows. 


As an example one considers the case when a external 


O' site of any polymer molecule reacts with a bound -CH^OH group 
of another polymer molecule. This results in five different 
- CH 2 “ linkages as shown below 

OH 


OH 


G H 


e-H 


u r 


h-i 


1 

on 


+ ai20H 


'T' 

O l! 


\ 


c> U 


+ CH20H 




Oj 


OH 


en 







+ H 20 


+ H20 


+ H20 


OH 

■^j^CH^( D>oH + H^O 


To be able to write the mass balance taking the reverse reactions 
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Into account It is necessary to distinguish the - CH 2 - linJcages. 
Therefore the following - CH 2 OH groups are distinguished, 

o VI 

(0 CH2OH CH2OH CH2OH 

(,-/N 

{ 
f 

All possible “CH 2 “ linlcages are shovm in table 5,1, With these 
the forward reactions beti-'jeen different species are shown in 
table 5,2, In table 5,3 all reverse reactions are given. It 
is assumed that represents the rate constant involving - CH 2 ~ 
and H 20 , 




/ c N' 


X 




A given bond can undergo several reactions for example 
1^2 linkage reacts with a water molecule in the following 
two ways ; 

Q) 


+ H2O + 20^ - 


+ H 2 O 


>0 + Og + p^ - p . 


Since it can react in two ways a rate constant of K^/2 is 
assigned to each. This vray the rate constant of each bond is 
decided. On the other hand when a - CH 2 " linkage react v/ith F, 
or the nature of the bond gets changed. For examx^ls a l^-l 

bond becomes 1,4 on reaction with F, or or ~ CH 2 OH, The 
foarward reactions involving the interconversion are shov/n in 
table 5,4. The state of different sites are not sho^irn in these 
reactions as it will lead to redundont counting. The reverse 
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reaction for inter con versions of bonds are sho^im in table 5 , 5 , 

One is nox'r in a position to write the mass balance 
equations for O*, 0 ., id., p , F, Q. , Q ' , p and balance equations 

Vb _L i X J, 

for all the lihl<:ages and five bound - CH2OH/ in the same xiray as 
done earlier. 

To get the exiaression for MWD, the balance equations 

for all the species P and Q are x^ritten. Terms due to reverse 

n n 

reactions are added to each of the equations derived earlier for 
the irrev-rsible case and the following equations results. 


d Lp 1 
""iftr™’" 


n -3 


_ CN 
SUM 

'-Vi* ( 

Djp + 

CQi !■ ) t i P 

‘-^l ^ ^_2 SUM n-J 


CM 

s^r* 

- LF2I &r 

1-2^ + 

(SK^^+k^) CpI L^n-l^ 


2 CN 

sum' 

cd 

CN 

SUM 

ilP^l ( X £ q ^\ + llQi! + 
2 

Coil ) 

-(n-1) 


K5 2Z- LQjlI Ch^o) +2K5 

CH20I 




i=n+l 




&.1, 

n ;^5 

(R-l) 


l=n 4 -l 



= (2K^+K4) Cp] ( [Qj^l + Top ) - -1^ tPj”' LfJ 



(R- ) 
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dt 


CM r „ 

SUM + 


- 9 ^ 

STJM 


LP2-i CQ|] + (2K^+K^) £p3 CQ23 


2 Of 
SUM 


^^3 3 - &3 ] ( + [Q^l + Ip^ ] ) 


mT =2 


■^^5 ^^3^ +Kg rH20] » [ q^] + 2 k ^ [H2OI Lp^l 


i=4 


1=4 


^5 t ^201 L Q3 ] 


(R-.3) 


d [;?! 

“dt sm r^sl tPlJ ^ ■*■ L?2 J ip23 

(2Ki^k,) 1: p3Cq 31 ^Igi^ CP 43 c f] 

“ “sum ^ ^ +£Qil ) - 3 K 3 £h 20 ]CpJ 

h-KsCh^oI 2J, [0^1 +2K5 (HzOl ZI CpJ + K3 roJCHjo] 

1 — b Xt -5 


It - = -( 2 I < i + K ,) fOnlOl - { on . CQ „ 1 + H ; L'Qnl ( ^ C (Q 

^ m=2 

+ CQ]1 +iP£l ) +LOn] [O'] +K^ CQil + t°i! )] 



CN 

$UM 


n-3 

. C Q^] C <SJv>-tl -'- 


1=1 


C M 

SUIM 


csiiot :® 


v\.ai 
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teac6H,,i.4<K^+K^) L'QiirQ^ ;i 


n-l' 


Cq {1 + 2 k^ cq 


1 ■' 




- K 5 fQ„\&20-i -K 3 (n_l) [H^o j + K 21, [Q “1 H,0 

i=n+l ^ 

n..:^/ 4 (R,^ 5 ) 

CqsI 

at = - ( 2 Kj+k 4 ) LQ 3 I Lp] - =n LQ 3 ] - si Psi ( ^C£ 4 l 

+ [Qll + [Oil) - COji (2Kj Cq'I + Kj^ CQj) ) 


** 2 Cl\F 

—m CP 3 JL 1 ''] + »^+K^) [DJCop + 2Kj^ rojl [0^3 

■*■ “sisr' LQ 23 ( [Oil + [Qil > - Kg - 2 K 5 P 33 

ALHjO] + K 2 Con CHjO] (r_ 6 ) 

A 


^ 2 Ji . . .,9 o 

= (K^+K^) i.Q^j + (3K^+K^)rQ^*]+2K^£Q«f 

+ “1^ rP2lfr] “ + 1^4^ £ Q 2 Hp 3 - CN (Q^l 

-Hr fts'l ^ S' ^'Qn? > 

m=2 

- [Q 2 I ( 2 K^ CQp +K^ CQ]} + ) - K 3 [■Q 2 \rH 20 J 


~K 5 G Q2irH20j + Kg Cq^I {H 2 O 3 (R-.7) 

jL*~3 
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Equation R-l is the rate exioression for n ^ 5* 

The reverse reaction in which disintegrates contributes the 

term (n-1) K^. Cp i !_H^O ] for there are (n-l) bonds in P , P 

can be produced .if any nM.) is attacfeed by a water molecule* 

Ther term Kj. /" ! C Q.l is thus accounted for. On the 

i=n+l 

other hand P^ is also produced when a P^ ( i ^ n+1) reacts with 

water. It can be noted that the breaking of bond from either side 

of P. (i!]^.n+l) will lead to formation of P and this leads to 

the term 2 k ^ C P.'j • In similar ways, the equations 

^ ^ i=n+l ^ 

for other species are written. 



TABLE 5«1 


Types of -CH 2 “L±ii]cages 


^ ^ ^ Q~Q lln Jcages 


OH OH 




OH OH 

I i 



OH OH 







(b) 0~p linlcages 
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TABLE 5.2 


FORWARD REACTIONS 


We are not calling a CH2OH as a. bond 
(a) R eact ions with P 

0 - 


4 K- 

1 . O* + F — j- 

e 


•> 


O' + O.+H^O 
e 1 2 


K- 


2 . 0 . 


3 , O^ + p 


IC 


K- 


© - Pe+Pi +H2O 

Q) 


Pi + p 


Pe + P 


+ H2O 


4 ) + H2O 


K 


■4 


© 


+ Oi- 0^+ H 2 O 


(b) 0 ^ with CH2OH 


© 


1. O ' + 

e 


2 , 0 ^ + 


3. 04 . (i) 


+ 0 .~ 0 '+Ho O 
1 e 2 

#• 0.-0 'tH^O 
i e 2 


4 . O^ + ^ 

5. 04 + (?) 


(c) 0 with CH^OH 


1 - V® 


2 , Og+ 



1 e Z 

-> (2;5) + 0.-0^+H20 

0 .- 0 *+H ^0 
1 e 2 


0 


b, (T© 


K, 


Oi |4 


C 

u>CH 2 

1 

^^CH 20 H 
CH2OH 

|)CH20H 

/ 

1 





3 . 


K. 



+ 
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0 *4“ 
e 

0 



0 4- 






5 /!) ■'•Pl-Pe+Ha® 


csM ; 

CH2OH <. 5 >< 3 Tt 

2,^ + Pj^-Pg+HjO CH^OH 


0. 

_A. 

with 

012 OH 

1. 

0, 

1 

+ ® 

2. 

0. 

1 

. © 

3, 

0. 

+ ( 3 ) 


1 


4. 

0. 

+ @ 


X 


5. 


+ ( 5 ) 


with 

CH^on 




K, 


K- 


K. 


K. 


^ 1 , 4 J + H2O 


,lv 5 ) + H2O 


^ 1 , 6 ) + H2O 


C)H 



an 




CH2OH 

e-Vi 

CH 20 H-f^ 


o H ^ 


(S) 

€ 5 ) 



i 

1 

i 

(y> 




©) 


+ ©) 


© 

- 4 - t /] ? 1 

K3 

C3) 



+ ( J) ■ 

,f 3 

Ci) 


CH2OH 




+ HjO 


+ H20 


+ H20 


+ H20 


+ H20 


+ H20 


k 

3 


7 
So 


m 

CH2OH Q 


H 

rz 


Ct' 

■4 


CH 20 H-^o> C'."' 
CH 20 H<]i^^"» 


0 V\ 


CH 2 QH 0 - 


.tl- CH^OH(^ 
2 




9- CH OH(t; 




■ 'h ' ^ 

'- 1 ^- CH 20 H-<^"’» 


Pe ■*■ 


K, 



+ 0.-0 +H0O 

1 e 2 


''■"10 


OH 

CH^OH {4 


K 


2 • Pe + lE) 
3 - Pe + @ 


4 


K. 


on 






4 - Pq + 14 


K, 


5. Pg + \5 


K. 


(272) + O^-Og+H^O -9 CH 2 OHQ 

{ 2,3 

Q 

(3/1 


+ 0 ^- 0 ^+H 20 

+ O .-O +H..0 *^0 CH.,OH-<£.>' •■’ 

JL ^ 


+ O.-O +H„0 
1 e 2 


tjl CH 20 H-<°^"'-' 


■ ^ 

i.j i. 

A iy-i.»> 
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(g) Rxn of P with p and CH-OH 


1. P + P 
P + P 


+ H 2 O 


Q» + H 2 O 


3. P +12) 


4, P +13 


5. P +( 


P +15 


^ + 0^+p^+ H 2 O 


^ + Oe-'Pe"- ^2° 




<2^ + (^q+Pq+ H 2 O 


+ 0g+ Pe+ f^2° 


7. P + 


8. 

P 

*Q 

A 

-> 

9. 

P 

+0) 








10. 

P 

+® 

4 

-*> 

11. 

P 

+(3) 

""4 



-> O-fSj + 20^ + H 2 O 


V<- / J- w y I V-/_ * XXr^\y 

V.1 — ^ e 2 

©) ^°i * «2° 


(h) Rxns of P X'7lth and 


1. P + Q| 


2. P + Q| 


3, P + Q, 


Pe + 20^ + 0^+ H2O + 


^ + 40| + H 2 O 




1/10) + 2p^ + 20^ + H.O 


4. P + Q, 


>2/12j + O^+p^+ 204 +H 20 


(i) Rxns of and with sites 


1- Oe Qi 


^ 0j_+20^-0^ + + H 2 O 
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(j) 


2. O- + 


3. 


0=. + Q-t 
e 1 


s. Pi + Qi 


6 . P^+ 


7. 0j_ + Q«^ 


8. Oi + Qi 


Pe °i 


10, p + Q 

1 


K, 


_fL 






Jll 




K3 









to 



> 1,9 + p^+ 0^+ O^- + H2O 

2,9 + 20^ + p^ - p^ + H2O 


+ Oq^ Pe+ Pi + H2O -p^ 

iSi) + =‘'S + H^o 


( 2 ,^ + 0„ + p + H„0 

\ 0 '*'0 2 

• O 1 » r-r 




K. 


K 


-> (^) + Og + Pg+ H2O 

3,4,; + O.-O^ + 20* + H^O 
vl..-'' 1 e e 2 


-“> Pe °i“ ^^2° 


Rxns of and Q| with CH 2 OH 



( 2 ^ + Qg + P^-Pi + ( 3 ) 

■i“> ( 1 ,^ + (T) + p^ 


K, 


K. 


2. + 2 


3, + 3 



2/lJ + C?} + Oj^ 

+ H2O - + Q 

1,2) + H20 + p_i_ + G) 


+ 0 + p 

e -^e 


i,i) + H20 + p^ + 

2,3) + H2O + + © 
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Qi + ( 4 ) 


K 


5. 


□l + 


6. QJ + (T 


■•■ ''3) 


8. + (s’ 

9. Q« + ( 4 ) 
LO. 


1 


{2 f6'j + H 2 O + Pj_ + 

■“> \3,2)+ H^o +0. + ( 2 ) 


IC 


^2 

^2 0 «• 0 ^ + + 0 ^ + p 


K 


2k, 


^ 2 ° + Pi + ( 1 ) 


2 ) + H 2 O + + ( 5 ) 


K 


'3 . r' 

t3,5_y + H^o +20 
2 k 

Tl '•5\ J. u 

2 


• - P. 

e -^1 

r ■'■•. 


. + (' 3^1 


X— V (i(^ + n o + 0. + (5) 


[^^ 2^ + H^O — 0 ^ 4 - 

2^1 ■" A 

' ‘,1/5J + H,0 + 0. +( 51 


11* Q^+ Q^ 


(i) 

2 K- 

-■■■- r- - ■■ iiA 

( 2 .!;) 

.In. 

to 

0 


■T* 

0 

2 K 

-.. .inTm,, 

(SA 

I 

to 

0 


"T 


- 2 - - 

'AA 

-i- 

H 20 


■' 

1 


2 K, 

A"A 





Cl > 9 ') 

-r 

H^O + 


A 



(A) + «2° + + Pi •" © 

1__7 ( 3 ^'^ + H 2 O + 20^ + 0j_ + ( 2 ) 


12. Q^+ Q^ 


13. Q^+ Q| 


2K 

2K 


^ “> fcz) + ^ 2 ^ • '^e' 

A* A 


+ 0e+ Ps+Pl+ A* 


4 K, 


i 

J 


■> ^. 7 ) + H 2 O + Og + p^+ 0 ^ +((2 

■:7 (2,1Q. + 20' + 0.+ 5 

VT' e 1 
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TABLE 5,3 
.REVERSE REACTIONS 


(a) With -CH2OH 


1 . 

2. 

3 . 

4 . 

5 . 


© 


K. 


+ H2O 


2 ) + H2O 



K 


0 

0 


+ H2O 


+ H20 


5 } + H^O 


K. 


K. 


F + 0 ^+p ~p. 
e -^e -^1 


20 * + F - O. 
e X 


0 ^ + F 


Pj_ + F 


Oe+Pe-V^ 


(b) With -CH2- Lirfcagej 



1,1) + H2O 


1^2) + H,pO 


K. 


0' + Oe + Pg-P. 




" o 


; + 204 - 0^ 


hKt 


+ H20 


Uf) 

1,4) -I- H2O 


Kr 


hi^c 


( 2 ) + 0^+ Pg-Pj. 

Qj * - °± 


'•“N 

1 / + 



C0 


+ H20 








Cl) + Oe+Pe-Pl 

2 ) + O. 


K. 


1 ,6) + H,0 


G) + 20^ - 

(ij * °i 
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^Z) + H 2 O 


hK, 


isK 


-7 Qi-Oe~Pe+Oe-*'Pe“I'i 

^ ^>P - 0^-Pg+ ( 1 ) 


J/8 3 + H 2 O 


(^1/ 9 ) '^^2^ 





"7 Qt- 0 ~p H- 0. 

1 e -^e 1 

hK^ 


3 

- Oe-Pe '+ P 

/7 



-> Qi-0e-Pe+20e-°i 




Kr 


■>>P 4 . 


© -0.- 


e Pe 


H 2 O V— 2q^ - 


\ 




hK, 




{2, 2 ) 


hK 


+ H 20 




- 7 { 1 ) + o^+p^- 0 . 

^ v_y e “^e 1 

"'■' G + Og+Pe-Pi 

^ -7^ + Oe+Pe-°l 






%K 


+ 20*-0. 
e 1 


^ /S] + o +p -o. 

t, ^ * £::i t 


“7 


e e i 


7 


(^2,4y + H 2 O 


GG + «2° 


^-> + 0 . 

h K 


y + Oi 


% K. 


-?-Q 


^ + Pi 


K 


“->© + 204 - 0 ^ 


ifIC 


-(i'- 


W + Pj 



57 


ni-i 


16 , 




17. f 


O ^ 

./v .V, 1'. 

(j ! < 






CH 



f;'H 



22, 


L ,4 


24. 


Oj-I 



O 


oy^i-i 

1 


\^2,6 ) + H^O 




■> © 


+ p. 

1 


^2 K 


'.2.7/ + H 2 O ^ 


Cl..] + Og+Pg-P^ 

^ Qi-Og-Pg+o^+p^-o^ 


'^2K, 


2,8 1 + H2O 


“> P -Oe-Pe-^ ^ ^ / 


+ P - 0 ^ - p^ 
e -^e 


2,9j + H 2 O 



© 

JSK5 

n . 


□i 



zj 

"V" 

Q-f 


■> P + 


© -20, 


C2.10)+ H 2 O (^ + p 


20- 




(2,11) 




+ H2O 


Q{ + 2o| - 0 204 

Qi „ 20* + 0 . 

X 0 1 


% K, 


-2 i^c 

. .. .. p «« 20 * 3 / 


hYi, 


(2,12j+ H2 


,0 


^ K, 


^ Q' - 20‘-0^-p„+P 
1 e e *^e 


Ql-20«-.0^-p^+P 



K 


3,1! + H.,0 


3,2 ) + H..0 
/ 2 



= ?©) 


^ + O +p -0. 

e 1 


+ Og+Pg-0^ 


%K5 


Hsj + p . 
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25. 

•: , _ J 

Gi':< ■' 

S Z v-^ 

) ^ 

V, .•-'' 

+ HjO 

^5 

->(4) 

26 . 

A. 

© 

+ H^O 

^^5 

“> Q{ 


e e -^e 






■ ■ ' '■)o i V 




28 • jl \Sy ■ '‘■^' 


.3^5 ) + MO 


o 

.3^6y + H 2 O 


P - 20| + 


d) 


'2K, 




K. 


Q‘ -20'+ p. 


-5? P - 20^ + k ) 


2Q' - 40' 

1 e 


O. 

X 
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TABLE 5.4 


Tn.'-t'srconvox'si.ons of Bo;id.i 


i,y + cx 


Stands for 


2 K 3 . 


sum of -CH20H,F,Q^,Q| 


''' ^2^ 


1,2 ) + oc 


(1,^ + H 2 O 


1,3) + CX 


(1,^ + CX 
Qi,5) + GC 
(1,7) + CX 


-> (Q) + H^O 


4-6) 


+ H 2 O 


■> (1,6 ) + H^O 

4ii} + «2° 


-> \l, 8 j + H2O 


l^-s) 


+ CX 


H 20 


(l, 4 j+ H2O 


1 , 9 ; + 


l, 5 j+ H2O 


^ (i-2) 


+ H^O 


l.ioj + CX 
• 


— > (l, 3 J + H2O 
( 1 ^ + H 2 O 


2 , 1 / + CX 


“> ( 1 / 9 ) + H2O 

->(^,6^ + H^O 


^2,3) + H„0 
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11 . 


f ^ 

K„ , 

V2^2j + GC 



12 . 

13. 

14. 

15. 


16, 


17. 


18, 


19. 

20. 


/' N 
2.3) + CX 




2,8) -I- C5C 


2 , 9 ') + GC 


(2,1^+ GC 


QC 



2 

'2^ 


K, 


K, 


K, 




K. 


K 


K. 


K, 

2K 

K, 

* 

2K. 

2k, 

K, 


K, 


2,4\ + H2O 



> ( 2 , 1 ) + H,,0 


-> \2,2) + H^O 


V2^8) + H^O 

\2,6j + tl^O 




■7 (2,5; + H 2 O 


“7 12,1) + H2O 


“7 \2,11/+ H^O 

^ / 

-7 C2/2; + H 2 O 


“7 ( 2 , 11 )+ H 2 O 



-> {2,3) + H2O 


-> \2^) + H2O 


■7 ^,9} + H2O 


-7 (j/l^^+ H2O 
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13,1' + CX 


-> ( 3,2 } + H «0 


2 \ + GC 
(”3,4'^ -h GC 


3,5^^ + GC 


(373) H2O 


.3,1) + H,0 


(3,5) + H2O 


(3,2j + H,0 


3,6] + GC 


{3,4] + H^O 
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table 5.5 


^(^erse Rxns in Interconve rsions 


1 . 

2, 

3. 

3. 

4. 


5. 


6 . 

7, 

8. 

9. 


(^,l) + H2O 

(1,^ + H 2 O 



+ H2O 
+ H2O 



(2)+ Hj< 


tl^7j+ HoO 


.l,8j+ H 2 O 


l,9l+ H^O 

V V ^ 




(1^ 


s 

! 
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10 . 


11 , 


12 . 


13 . 


14 . 

15 . 

16 . 

17 . 

18 . 
19 , 



+ H2O 







+ H2O 




f 
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aiAP':ffiR 6 


CONCLUSIO.WS 

A kinetic model for irreversible Wovolac type phenol 
formaldehyde polymerization has been presented and equations 
governing the molecular v;eight distribution for batcli as well as: 
hanogeneous continuous - flow reactors have been derived, - The 
sot of non-linear differential equations for batch reactor are 
solved by Runge-Kutta method of order 4, Brovm's algorithm 
is used to solve the sot of non-linear algebraic equations : 

cierived for HCSTRs, This is found to be more efficient than 
the Newton’s technique of solving non linear algebraic equations, i 
The Mv’JD thus obtained fior- batch reactors andi HCSTRs as a function : 
of time ^ro presented, 

R1,R2,r 3 and Cp 3 ^ / Lp 3 ^ axe four parameters, vhosQ i 
effect upon the I'lWD of the polymer have been examined. The M^JD 
were found to be relatively insensitive to parameters R1 and R2 
and vKiie extremely sensitive to R3, For a given R3, the I 

distributions of as x^ll as fall monotonically and the [ 

conv-.-rsion of phenol in HCSTRs is smaller and □it ! 

larger than those for batch reactors. For a given residence tirae, 
the number average molecular weight, for HCSTR is lower but | 

but the polydispersity index, -f' , is larger. The variation of v i 
R3 has jess proEiounced effect upon the liWD from HCSTR than that 
from the batch reactors. 
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